load my_signal_with_noise5.mat;

% EMD 分解和去噪  
imfs = emd(signal_noisy); % 经验模态分解
num_imfs = size(imfs, 2);

% 选择低频分量（前几个IMFs）重构信号
y_denoised = sum(imfs(:, 1:2), 2); % 只取前两个IMF分量作为重构

% 绘图
figure;
subplot(3,1,1);
plot(t, signal,'b');
hold on;
plot(t, y_denoised,'r');
hold off;
title('蓝色为原始信号，红色为处理后的信号');
xlabel('时间 (s)'); ylabel('幅度');

subplot(3,1,2);
plot(t, signal_noisy);
title('带噪声信号');
xlabel('时间 (s)'); ylabel('幅度');

subplot(3,1,3);
plot(t, y_denoised, 'r');
title('EMD去噪后的信号');
xlabel('时间 (s)'); ylabel('幅度');
